Longitudinal and spin/ valley Hall optical conductivity in single layer M0S2 
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A monolayer of M0S2 has a non-centrosymmetric crystal structure, with spin polarized bands. 
It is a two valley semiconductor with direct gap falling in the visible range of the electromagnetic 
spectrum. Its optical properties are of particular interest in relation to valleytronics and possible 
device applications. We study the longitudinal and the transverse Hall dynamical conductivity which 
is decomposed into charge, spin and valley contributions. Circular polarized light associated with 
each of the two valleys separately is considered and results are filtered according to spin polarization. 
Temperature can greatly change the spin admixture seen in the frequency window where they are 
not closely in balance. 

PACS numbers: 72.20.-i, 75.70.Tj, 78.67.-n 



I. INTRODUCTION 

Since the isolation of graphene^ the search for new 
two dimensional atomic membranes possibly with novel 
functionalities has intensified. M0S2 3 in its single lay- 
ered form is a two valley direct band gap semiconduc- 
tor with gap in the visible and so is of interest for de- 
vice applications. Through excitation by right and left 
hand polarized light excess populations of a selected val- 
ley can be generated which make this material ideal for 
valleytronics Because the two inequivalent valleys are 
separated in the Brillouin zone by a large momentum, 
intervalley scattering should be small. Consequently the 
valley index becomes a new degree of freedom analogous 
to spin in semiconductors. Just like manipulating spin 
has lead to spintronicsj^ manipulating valley index can 
produce new effects including using it to carry informa- 
tion. As an example, in the context of graphene Xiao 
et.alj 1 ^ showed that a contrasting intrinsic magnetic mo- 
ment and Berry curvature can be associated with carrier 
valley index. As a second example a valley filter device 
is described by Rycerz et.al.— 

An important issue is possible pathways to achieve val- 
ley polarization i.e. populating states preferentially in 
one valley* 1 ^ One way discussed theoretically ^ 10 ' 13 ' 14 is 
to use circular polarized light and this has been demon- 
strated recently by three experimental groups^— ^ Se- 
lection rules on the absorption of right (or left) handed 
polarized light assures that this radiation excites almost 
exclusively charge carriers residing in a single valley with 
index t = +1 (or —1). The physical quantity that comes 
into the description of such processes is the real part of 
the AC optical conductivity for right (+1) or left (-1) 
polarization Rea± as a function of photon energy lj. In 
terms of longitudinal <j xx (oj) and transverse o~ xy {uj) con- 
ductivity a±(ui) = <j xx {uj) ± i<T xy (uj). The conductivity 
also plays an important role in determining optical prop- 
erties of nanostructures as we can see from the references 
|l5Ml6lll7lll8Hl9ll20l ]. In order to solve the Maxwell equa- 
tions for systems with graphene sheet between two me- 
dia with different dielectric constants, one has to know 



the conductivity of two dimensional graphene. Results 
for the conductivity of single layer M0S2 will be use- 
ful if similar devices were built from M0S2 instead of 
graphene. 

M0S2 is a layer of molybdenum atoms between two 
layers of sulfur in a trigonal prismatic arrangement which 
does not have inversion symmetry. In momentum space 
at the K and —K points of the honeycomb lattic o 21 ' 22 the 
valence and conduction bands are separated by a large 
semiconductor direct band gap and there is a large 
spin orbit coupling leading to a spin polarization of the 
valence band, spin up f and down J, as the z -component 
of the spin operator s z commutes with the Hamiltonian 
and hence remains a good quantum number. A mini- 
mal Hamiltonian which describes the band structure of 
M0S2 (valid near the main absorption edge)is found in 
reference 4 with parameters based on first principle cal- 
culations for the group- VI dichalcogenides i 2 1 ' 24 A more 
complete theory is found in the reference [l4[ where a 
calculation of the dc Hall conductivity and Berry curva- 
ture over the entire Brillouin zone is presented. Part of 
the Hamiltonian describes the dynamics of massive Dirac 
fcrmions which are known from the graphene literature 
to have an optical response quite different from that of an 
ordinary 2D electron gas. For example a universal con- 
stant background conductivity of e 2 /(4fc) 25 i 26 (for a single 
spin and single valley this constant becomes e 2 /(16h)) 
is predicted and observed for photon energy oj greater 
than twice the chemical potential /j, for massless Dirac 
fermions. We use this Hamiltonian to calculate the dy- 
namic optical conductivity as a function of photon energy 
to several electron volts. Results are presented for longi- 
tudinal <7 xx (lo) as well as transverse conductivity a xy (uj) 
which is separated in charge, spin and valley Hall con- 
ductivity. Appropriate contributions from different spin 
channel are presented separately and special attention is 
payed to the effects of right and left handed light polariza- 
tion. Temperature effects are also considered. In section 
II we present the Hamiltonian and the Green's function 
on which our calculations are based. In section III the 
mathematical expressions for the conductivity are evalu- 
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ated and results for spin and valley Hall conductivity are 
given for several representative values of the chemical po- 
tential. Results for circular polarized light are found in 
section IV and section V contains a summary. 



II. FORMALISM 

The Hamiltonian for M0S2 at K and —K points is 

1 



Hq = at(rk x a x + k y a y ) + y<r 2 - ^-^Y 



(1) 



with 2 A the spin orbit splitting at the top of the valence 
band and we take 2A = O.lbeV, a the lattice parameter 
3.193^4, t the hopping t = l.leV, a the Pauli matri- 
ces and S z the spin matrix for the z-component of spin 
s z which is a good quantum number here. The index 
t = ±1 is the valley K(—K) respectively and A is the 
direct band gap equal to 1.66ey between valence and con- 
duction bands. The x and y velocity components based 
on |T]) are 



8Hq at 
hdk x 



-T(J X ,V V 



8Hq at 

hdk v h 



(2) 



Here we will be interested in charge, spin and valley 
I 



current given respectively as j = ev, j s = |s z v and j" 
tv. The eigen energies and vectors of ( [1]) arc 



and 



u n (k) 
atk x 



(3) 



-Ay2±ya2t2 fc 2 + (A72)2 rkx+tky 
at P 



a 2 t 2 k 2 + (-A'/2 ± y/a 2 t 2 k 2 + (A'/2) 2 ) 2 



(4) 



with A' = A - Xts 



For later reference the Berry 



curvature for the conduction band is 



O c (fc) = z ■ V fe x (u n (k)\i\7 k \u n (k)) 
2a 2 t 2 A' 



[A' 2 + 4a 2 t 2 k 2 } 3 / 2 



(5) 



with valence band ft v (k) = — O c (fc). With these solutions 
the green's function Go(k, iui n ) with Matsubara frequen- 
cies oj„ is 



(iftwn + n - \ts z /2)I + (—)a z + at(rk x a x + k y a v ) 



(ihw n +fi- Xts z /2) 2 - a 2 t 2 k 2 - (Al)2 



Gj(k, iw n )I + G z a z + G x a x + G y a y (6) 



which gives 



Gj(k, iw v 
1 



2 ihion +fj,- Xts z /2 - Ja 2 t 2 k 2 + {^-) 2 



where is a sum over momentum, Im means taking the 
imaginary part. The longitudinal conductivity ct xx (lu), 
charge Hall conductivity <j xy {oj), spin Hall conductiv- 
ity a xy (ui) and valley Hall conductivity a xy (u>) are given 



ifkv n +n- Xts z /2 + Ja 2 t 2 k 2 + (Al)2 



G z = 
G x = 
G v = 



A' 
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(ihw n ■ 


f H- Xts z /2) 




atrk x 




Jrfl- Xts z /2) 




atk y 


(ihw n - 


f n- Xts z /2) 



G I (k,iu n ) 
Gi(k,uj n ) 
Gi(k,iu n ). 



The density of state N(e) is given by 

N(e) = -— i"mG 7 (k,io; n - > e + id). 

7T — ^ 



(7) 



(8) 
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iuh 2 ^— ' 

k,i 



Tr(a x G (k,iuJi)a x G (k,iuJi + iu n )) 



i(jj n — >uj-\-i5 



cr xy (ui) 



e 2 a 2 t 2 



-T 



iujh 2 ^— ' 
k.i 



Tr(Ta x G (k,icJi)a y G (k,iwi + icJ n ))iu n ->u+iS 



k,i 



Tr(-s z Ta x G (k,iuJi)a y G (k,iuji + iu) n ))iu n ->u+i6 
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<J xv\ UJ ) = ""TT^ / 

Tr(a x G (k,iu>i)a y Go(k,iLOi + «'w n ))iw„->w+« 
After simplification we get 

„2 



(9) 



0-1/ H = E - /(£-)]<? (10a) 



with 



(4a 2 i 2 fc 2 + A' 2 ) 



(/iw + i<5) 2 - (4a 2 t 2 A; 2 + A' 2 ) 



(10b) 



for spin and valley Hall conductivity, g in (|10b[) is to be 
replaced by s z g/e and rg/e respectively. Here f(x) is 
the Fermi Dirac function, which contains the chemical 
potential /i. The longitudinal conductivity is 



2a 2 t 2 (4:a 2 t 2 ky + A' 



h 2 ^ iuj\A' 2 
[/(£+) -/(£?")] 



. 4a 2 t 2 fc 211/2 



(ftw + ?5) 2 - (4a 2 i 2 fc 2 + A' 2 ) 



(10c) 



To obtain (|10c[) we have included only the interband 
transitions. There is an additional intraband contri- 
bution which provides a delta function contribution at 
uj = 0. When residual scattering is included in the calcu- 
lations, the intraband piece broadens into a Drude peak 
which can overlap with the interband contribution. But 
in the pure limit which is the case we are considering here 
we need not consider this contribution. What replaces 
the Berry curvature in the expression for the longitudi- 
nal conductivity factor 



2a 2 t 2 (4a 2 < 2 fc 2 + A' 2 ) 
H " ~ iuj[A> 2 + 4aH 2 k 2 } 3 / 2 



(lOd) 



For k — this factor in a xx and a xy agree and as we will 
see later this leads to a valley selection rule for circular 
polarized light; at finite k however the cancelation is no 
longer exact. 



III. RESULTS FOR SPIN AND VALLEY HALL 
CONDUCTIVITY 

In Figured] we show our results for the Hall conductiv- 
ity cr X y(tL)) vs. hu in units of e/h for the four values of the 
chemical potential shown in the inset of the top frame 
where the band energies E^(t,s z ) are sketched (equa- 
tion ([3])). n = — 0.995eF falls below the top of the lowest 
energy spin polarized valence band; \i = — 0.845eV" falls 
between the two valence band; /i = 0.09eV is between va- 
lence and conduction band(insulator) while fi = 0.92eV 
falls in the conduction band also spin split but this split- 
ting is small. In the main frame which has 4 frames each 




1 2 3 

ftco(eV) 



FIG. 1. Real (solid) part of the Hall conductivity a xy vs. u, 
and its imaginary part (dashed) for rs z = 1. The dotted and 
dashed dotted curves are for rs z = — 1 instead. There are 4 
frames from top to bottom for four values of chemical poten- 
tial fx namely -0.995, -0.845, 0.09 and 0.92eF respectively. 
The inset provides a sketch of the bands in M0S2 and how 
the four values of fi relate to these. In all cases the bands are 
cut off at momentum ka = 3.0. 
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FIG. 2. The energy bands (solid) for M0S2 with chemi- 
cal potential jj, level indicated by a dashed horizontal line. 
The extremal optical transitions are also indicated in each 
case. From top to bottom frame fj, = —0.995, —0.845, 0.09 
and 0.92eV. 



for a different value of fj, we show separately the con- 
tribution for ts z = ±1 with a xy (u)) the sum over both 
index s z (spin) and r (valley). The solid line is the real 
part for ts z = 1 while for ts z = —1 the dotted line 
applies. The dashed and dash-dotted are for their re- 
spective imaginary part. We note a sharp onset in these 
two last curves the energy of which can be traced to the 
minimum energy associated with possible interband op- 
tical transition as shown in Figure [2j As an example, in 
the top frame we see that, for /i = — 0.995eV, the onset 
of the interband transitions for spin f occurs at higher 
energies than for spin J, (here the valley index has been 
taken to be 1 ). Corresponding to the onset in Ima xy 
there is a peak in its real part at this same energy, as 
they are related by the Kramers-Kronig relations. The 
results presented were obtained through numerical eval- 
uation of equation (|10a[) . The numerical values do have 
some dependence on cut off used on the energy (atk) 2 . 
Here we have set the cut off ak to be 3 and restrict our- 
selves to photon energies below 3eV. In this energy range 
choosing a larger cut off makes no difference to the re- 
sults which have converged. For photon energies above 
3eV, a range not considered explicitly in our figures, in- 
creasing the cut off even to infinity has no qualitative 
effect on the results for the imaginary part of a xy {w). 
It affects the real part more. With a cut off, there is a 
zero in Rea xy (uj) at some high energy. As the cut off is 
increased, the energy of the zero in Rea xy (uj) moves to 
higher energies and in the limit of infinite cut off no zero 
remains for finite to as it has moved to infinity. Taking 
the infinite band limit has the advantage that simple an- 
alytic expressions can be obtained which can be useful. 
For example it is straight forward to show that 



Ima xy (uj) 



2(A — A) 



16h fkj 
A haj A frbj 

x [/( 77 + ir ~ m) - /(o - ir - 0) 



9(hiu — (A - A)) 



(11) 



where we have made explicit the value of chemical poten- 
tial attached to the thermal function f(x) (Fermi Dirac 
distribution f(x) = l/(e l3x + 1) with f3 = l/T). Note the 
onset at (A — A) and the 1/lu drop in this function as a 
function of uj. While there are some quantitative differ- 
ences of the form (|f Ij) with our numerical results there 
are no qualitative changes. Another reason for restrict- 
ing the range of photon energies considered, as we have 
done here, is that the model Hamiltonian (I) is itself valid 
only near the main absorption edge. A first principle cal- 
culation of the dc Hall conductivity and Berry curvature 
over the entire Brillouin zone which goes beyond what we 
have done is given by Feng et.alJ^ Here we really cannot 
access accurately the high energy region. 

To get the spin and valley Hall conductivity from the 
results presented in Figure [T] components need to be 
added according to the weighting hs z /(2e) and r/e re- 
spectively as noted in relation to equation (|10a)l . Results 
are presented in Figure [3] The real part of the valley 
Hall conductivity is the solid curve with dashed the cor- 



5 




0.6 
0.3 



0.0 



04 


-0.3 



O 0.6 



0.3- 



-0.3 
0.6 

0.3 

0.0 h 

-0.3 



T=10K 



H=-0.845eV 




o.o — ■-- 



H=0.09eV 




H=0.92eV 



1.2 1.5 1.8 2.1 

#w(eV) 



2.4 



FIG. 3. The real (solid) and imaginary (dashed) part of the 
valley Hall conductivity a xy vs. ui. The dotted and dashed 
dotted curves are similar results for the spin Hall conductiv- 
ity case. From top to bottom frame the chemical potential 
H = -0.995, -0.845, 0.09 and 0.92eV. The inset in Figure 
[T] provides a sketch of the bands in M0S2 and how the four 
values of \x relate to these. 



responding imaginary part, while the dotted and dashed 
dotted are for the spin Hall conductivity. The differences 
in weighting s z and r has a large effect on the resulting 
Hall conductivity. For example the real part of the val- 
ley Hall conductivity is everywhere positive in our model 
while the spin Hall conductivity starts negative at lu = 
and is very small in comparison to the valley Hall conduc- 
tivity. Both show peaks at the onset energies associated 
with their imaginary part. Because of the spin splitting 
of our bands there are two peaks seen clearly in the real 
part of the valley Hall conductivity (solid curve), after 
which it drops off gradually with increasing u). Conse- 
quently the real part of the spin Hall conductivity (dotted 
curve) changes sign with increasing w, because the two 
peaks found in Figure 1 (solid and dotted curves) are as- 
sociated with opposite spin, this sign changing behavior 
has also been found in the Figure 4(a) of the reference 
[14] . While the results show quantitative variations with 
value chemical potential (see inset in Figured]), there is 
no qualitative change. 

The zero energy limit (u> = 0) of the Hall conductivity 
is of interest and can be worked out analytically in the 
infinite band limit. At zero temperature (T = 0) the 



results are (a xy m un it of |, <J v xy in unit of | 



Kyi" 

O = 0) = 



2/i 



7T A 2 - 4/i 2 

1 2//A - A 2 



7T A 2 — 4/i 2 

for ji in the lowest valence band; 



1 A - 2A + 2/i 



and 



1 A -2/i 



2/i 



(12a) 
(12b) 

(12c) 

(12d) 

7T 4/i 2 - A 2 7T 4/* 2 - A 2 ^ 12e) 

for /i above the conduction band. Note in particular for 
/i between valence and conduction band T xy (uj = 0) = 
and o v xy {uj = 0) = Here we do not have a spin Hall 
insulating state. Note that for the real part of the Hall 
conductivity at zero temperature an analytical expression 
exists, for example for r = 1, s z = 1, 

„2 



2tt A - 2/i 2?r A 

for /x between f, \. valence band; 

and — 

for /1 between valence and conduction band and 

A A — 2/i ,1 2uA — A 2 
and 



Rea 



x In 



hiU — 2-y/Xmin hid 



(A -A) 
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with 



max[(/i 



A 2 

2> ' l 2 



(ka) 2 cut t 2 + ( 



A-A, 



(13a) 

(13b) 
(13c) 
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(ka) cu t is the cut off for ka, in the infinite band ap- 
proximation {ka) cut — > oo while in our numerical results 
(ka) cu t = 3.0. Similar results can be found, for example, 
in the reference I31L 



IV. CIRCULAR POLARIZED LIGHT AND 
TEMPERATURE EFFECT 

For circular polarized light the appropriate optical con- 
ductivity is 



a±(u) = a xx (uj) ± ia xy (uj) 



(14) 



where the longitudinal and charge Hall conductivity 
given by (|10a|10c[) . Results for the absorptive part of the 
conductivity Rea±(uj) are presented in Figure 2] where 
we show separately Rea xx (solid curve) for s z = +1, 
and dotted curve for s z = — 1, to be compared with the 
dashed curve for Ima xy with s z — +1 and dash-dotted 
for s z = — 1. For a given spin the onset in each pair 
of curves for longitudinal and Hall conductivity are the 
same. The two begin to deviate from each other as to is 
increased above the threshold where the Hall conductiv- 
ity falls below its longitudinal value. 

For the infinite band case we have seen in (fTTj) that 
Ima xy drops like 1/uj. Similar algebra for the infinite 
band limit gives for the longitudinal conductivity 



Rea xx (uj) = - 
, A Hui 

x[/(x + ^r 



e 2 



(A-A) : 



k 2 UJ 2 

, , A 

^2 ~ ~~2 



]6(Huj 



(A -A)) 



(15) 



which drops less rapidly with increasing uj than does 
Ima xy (uj) in equation (|lip and provides a check in our 
numerical results shown in Figure 2] for a finite cut off 
ka = 3. 

The results in the top frame of Figure |4] apply to 
temperature T = 10K, while those in the bottom are 
for T — 300X. For the specific value of /i chosen i.e. 
fi = — 0.845eV the chemical potential falls between the 
two spin polarized valence bands. On comparing with 
the top frame we note considerable temperature smear- 
ing of the spin up band. But the spin down band by 
comparison is much less affected as it does not fall at the 
Fermi energy but is everywhere below fi and hence does 
not respond to temperature as effectively. 

In Figure[5]we present our results for Rea±(u>) vs. uj in 
units of e 2 /h for valley r = — 1 in the case /i = — 0.845ey. 
The top frame applies to T — 10K while temperature is 
300K in the bottom frame. Results are presented sep- 
arately for s z — ±1. We note that Re<j + (uj) is large 
and comparable in size for cither spin up or down while 
Re<j-(uj) is very small in comparison. This is under- 
stood from the optical selection rule which apply to a xx 
and a xy . In the infinite band limit we can show that 
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FIG. 4. The absorptive part of the longitudinal {Reo xx (u))) 
and transverse Hall conductivity (Ima X y(uj)) as a function of 
photon energy u). The chemical potential is n = —0.84561/. 
The top frame is for temperature T = 10K and the bottom 
for T = 300K. The solid and dashed curves are for Re<r xx {uj) 
and Ima xy (u}) with r = +1 (valley index) and s z = +1 (spin 
index), while dotted and dashed dotted are for s z — — 1. 



Rea xx (oj) =F Ima xy (uj) 



d(huj - (A - A)) x 



[1± 



( A - A )i2r,, A 



f(-, 



hiu . , A hiu . , . , 
T -U)~f{^ - T -M)](16) 



which implies a perfect cancelation at the onset energy 
uj = (A — A). Here the equation (fTo) is for s z = +1 
band with —A — > A for s z — —1 case. Effectively light 
polarization provides valley selection to a very good ap- 
proximation. It is exact at the onset energy and remains 
quite good even at uj = 4eV, for which energy the light 
polarization will reduce the absorption from this valley 
to about 20% its value for the opposite polarization. The 
polarization selection rule found here for the dynamical 
conductivity agree fully with those previous discussed by 
Yao et.ali^ for inversion symmetry breaking Hamiltonian 
in two valley systems, where they relate the optical selec- 
tion rule to the orbital magnetic moment and the Berry 
curvature of particular valley considered. It is interesting 
to look at the spin up spin down admixture of the optical 
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FIG. 5. Results for the absorptive part of the conductivity for 
circularly polarized light, g±{uj) = a xx ±ia xy . Here the valley 
index r = — 1. The solid is for Rea+ the dashed for Rea- 
both for spin index s z — +1. Fors z = —1 the dotted curve 
is a+ and the dashed dotted for <r_. The dominant response 
in this valley is from <r+ right handed circular polarization, 
also shown on the same figure is the up to up plus down spin 
polarization for a+ (solid red curve). The top frame is at 
T = 1QK and the bottom at T = 300^. 

response Rea + (uj) vs. w. We define 

spin up/ \ 

P( W ) = — ^± &l (17) 

v ' spm up / \ . spzn down / \ v / 

The results are shown as the red curve in Figure [SJ We 
note that for most frequencies P(w) is close to 1/2 but 
that around the onset energy for er + there is a region 
where P(w) falls below this value and can be close to 0. 
This is easily understood as a direct consequence of the 
displacement in onset between Pecr+(oj) for spin down 
(dotted curve) and for spin up (solid curve) which implies 
a deficit of spin up electron in the region between these 
two onsets. Temperature can have a large effect on the 
position and shape of P(w) in the region of the onset as 
seen in the lower frame of Figure [5j This can be traced 
to the fact that the spin up electrons are much more 
susceptible to temperature smearing for the case consid- 
ered here as we have already noticed. The solid curve 
(t) for 0+ now extend to lower energies than does the 
short dashed curve (4-). Temperature can in fact change 
the magnitude of P(u>) in the region of interest from less 



than 1 /2 to larger than 1/2 with this entire region shifted 
toward lower energies. The peak in P(w) for T = 300K 
is also broaden as compared with the valley in this same 
quantity for T — 10K. The spin admixture in this region 
can be manipulated with temperature. 



V. SUMMARY AND CONCLUSIONS 



We presented expressions for the dynamic conductiv- 
ity of M0S2 based on a simplified Hamiltonian which 
includes spin orbit coupling and band structure param- 
eters fit to first principle calculations. Using the Kubo 
formula, the final expressions reduce to a sum over mo- 
mentum k centred about the two valley points K and 
—K defining the corners of the honeycomb lattice in the 
Brillouin zone. The bands are spin polarized and a sum 
over spin and valley appears. The transverse conduc- 
tivity is split into charge, spin and valley Hall conduc- 
tivity and, in all cases, depends on an overlap of the 
Berry curvature multiplied with a common sum of two 
energy denominators linear in frequency plus appropriate 
temperature factors and channel dependent indices. The 
expression for the longitudinal conductivity is not quite 
as simple but reduces to the common form around k = 
(zero momentum). Simplified analytic expressions can be 
obtained in the infinite band limit. We have presented 
numerical results for all these quantities as a function of 
photon energy separating out, in each case, the contribu- 
tion from the separate spin channels. In the numerical 
calculations we use a cut off on momentum of ka = 3.0. 
The results are not qualitatively different from the infi- 
nite band case although there are important quantitative 
differences. The effect of temperature is considered. It 
provides a smearing that tends to obscure the separate 
contributions to the conductivity. It also can change very 
significantly the spin admixture in the frequency window 
just above the main absorption threshold where there 
can be an important imbalance between up and down 
spin, for a given valley when circular polarized light is 
employed. The absorptive part of the conductivity for 
right and left handed polarized light, which is the appro- 
priate quantity for valley tronics is also computed, and it 
is found that second valley contributes very little to the 
absorption. For example it is less than 20% at 10 = 4eV. 
with main absorption peak at 1.7eV. The DC limit of 
the Hall conductivities is obtained analytically. For the 
parameters appropriate to M0S2, the real part of the 
valley Hall conductivity is positive and large in value as 
compared with its spin Hall counterpart which is nega- 
tive for some values of the chemical potential. We hope 
that our calculations can help further our understanding 
of the optical properties of MoS^- 
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